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Dark soliton (DS) and Josephson vortex (JV) in quasi-lD long Bose Josephson junction (BJJ) 
can be interconverted by tuning Josephson coupling. Rates of the interconversion as well as of the 
thermally activated phase-slip effect, resulting in the JV switching its vorticity, have been evaluated. 
The role of quantum phase-slip in creating superposition of JVs with opposite vorticities as a qubit is 
discussed as well. Utilization of the JV for controlled and coherent transfer of atomic Bose-Einstein 
condensate (BEC) is suggested. 
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Solitons, in general, and the DS in particular, 

continue to be a fascinating subject for study. In ID the 
snake instabihty 0] is suppressed and the DS is a stable 
particle-like object (apart from the slow phonon induced 
decay Q)- Vortices in atomic BEC have been studied 
in great details theoreticaUy and experimentally as well 
(see in ref.Qj). 

Superfluid current circulation can exist in two paral- 
lel quasi- ID waveguides, coupled by a uniform Joseph- 
son tunneling 7 > along their length (see Fig.l), akin 
to the JV in superconducting long Josephson junction 
Traditionally 0, phase variation only is considered 
within the frame of the Sine-Gordon (SG) equation. As 
it turns out, such description is insufficient for the BEC 
waveguides in the quasi-lT) regime, where the ID Gross- 
Pitaevskii (GP) equation in the axial direction should 
be employed. Short BJJ have been thoroughly exam- 
ined in refs.0, H, 0- Two coupled waveguides were al- 
ready studied as well with focus on bright solitons [lo|. 
Quasi-ID BEC were created in a variety of magnetic 
traps 01 . Intrigu ing perspectives are offered by BECs 
on microchips |l2t ll,*^ . Coherence between two parallel 
elongated BECs has been demonstrated in the seminal 
MIT experiment Il4l and is currently being under in- 
tense investigation [iSlll^ . Two parallel waveguides with 
varying separation from each other were designed as an 
interferometer in refs.0,0|. As noted in ref.^3|i the 
considered s etup has precise optical analogy - dual-core 
optical fiber [l7| . Similar system with point-like coupling 
has been considered in ref . [13 with emphasis on possible 
application in atomic interferometry. 

Here we show that as the ratio v — j/fJ., with /i be- 
ing the chemical potential, becomes smaller than some 
critical value the DS transforms into the JV sponta- 
neously. This breaks the time-reversal symmetry. Con- 
versely, as v exceeds I'c, the JV transforms into the 
DS, which restores the symmetry. The DS^JV inter- 
conversion effect is a reversible ID analog of the 3D DS 
snake instability 0|- In contrast to the 3D, where the 
DS irrecoverably decays into vortex rings, the DS in the 



quasi- ID BJJ can be controUably restored from the JV 
by tuning above the critical value. The thermal and 
quantum phase-slip effects can restore the symmetry as 
well. 



The initially closed vortex circulation 
opens when the vortex encounters a 
boundary. The redirected BEC flux 
transfers BEC atoms from Reservoir 1 
to Reservoir 2. 
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FIG. 1: Schematic representation of the long BJJ. 



Model: Action. Dynamics of the bosonic fields Vfej 
/c = 1, 2, is controlled by the action 



S— dtdxC, £ = £1 + £2 + ^^12, 



(1) 



where t is time and x is the axial coordinate; the La- 
grangian density C consists of the parts 

Ck = ^ihi^lij,] - l^lV^fcp + /ifel^fcP - ^\^k\\ (2) 
2m 2 

describing the dynamics along well as the Josephson 
tunneling 



7'0iV'2 + c.c; 



(3) 
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chemical potentials /Zfc = /z in each waveguide will be 
kept identical; m is atomic mass; g = Anh^ a / [mr^) — 
the effective ID interaction constant with a > be- 
ing 3D scattering length and r_L standing for the effective 
width of the waveguides. Strictly speaking, the standard 
Josephson coupling Q is valid for <C 1. As ^ 1, 
higher order terms in 'ipi^2 should be considered. How- 
ever, while such terms may change the numerical value 
of Vc, they will not eliminate the interconversion effect. 
Thus, in what follows, we will not discuss the higher or- 
der couplings. 

Model: Dissipative function. Significant damping of 
Josephson dynamics can occur due to tunneling of the 
normal component ^] between two BECs. Phenomeno- 
logically, the dissipation can be introduced through the 
dissipative function Fd, so that 



d SC 
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in accordance with the standard procedure |20|. In gen- 
eral, Fjj must be positively defined function of the time 
derivative of physically observable quantities We 
choose it in the minimal form 



Fn = / dx^, 



(5) 



the total number of atoms. We note, however, that they 
violate the Galilean invariance (given by the transforma- 
tion dt dt — ipk exp[i(yx -I- V'^t/2)]ipk, with 
V being the velocity of a new frame moving along x). 
In this regard, it is important to realize a limited na- 
ture of the above phenomenological approach — it can 
only be applied in the case when the tunneling of the 
normal component does not conserve linear momentum 
(along the waveguides), that is, when scattering on im- 
perfections of the trapping potential is significant as in 
ref.lij. 

Dark soliton and Josephson vortex. Static one soliton 
solutions of eqs. H6I7() in infinite medium belong to a 
family 



1,2 



\/l + ly tanh(px) ± i 



cosh(pa;) 



(8) 



The DS, characterized by = -02 corresponds to 

p = a/T+I^ and B = (that is, total density has zero). 
The static JV satisfies combined symmetry — time re- 
versal and reflection {ipi — ijj, ■02 = and, thus, the 
equation 



(9) 



can be taken 



where p = IV'iP — IV'2p- AH the information about 
normal component is included into the kinetic coefficient 
a, which can be related to the dissipation rate of small 
Josephson oscillations considered in ref. i^. Employing 
eas. Hll5|) for uniform p and small relative phase, we ob- 
tain damped Josephson oscillations p + uP'jp + K/j = 0, 
where the Josephson frequency and the damping coef- 
ficient are = 2yJ^{p, -f 2^)/h and k = 87(^ -I- 
j)/{h?ga), respectively. The value of k, which deter- 
mines a typical relaxation time t ^ a, 
from the microscopic analysis 

We introduce the units of length 1^ = h/^ 
relation length), with no — p/g being the average ID 
density in a single uncoupled (7 = 0) waveguide, and 
of time to — h/p. Then, setting ipk — > ^/noipk, the 
action (^) becomes S = HSq J dtdxL, where 5*0 = Icn^ 
(the validity of the GP regime is justified by 6*0 S> 1 
|3); L = Li + L2 + ii2, with Lk = i(0feV'fc + c.c.)/2 - 
|VV'fc|V2 + l^fcP - 1^-^172 and L12 = j^V'i*V'2 + cc. The 
dissipative function density takes the form hS^fP' /2a, 
with a = hat^/nQ. Employing eas. (|ll5| ) in these units, 
we obtain 

1 

#1 - -/O^l = - V'l + l^lPV'l - t^V'2; (6) 

a Z 



ii'2 + —P'^2 

a 



-■4'2 - tp2 + \tp2\'^'lp2 - l^^l- (7) 



The dissipative terms ~ p resemble the phenomeno- 
logical dissipation introduced in ref. '21*1. These conserve 



Its solution © is given hy p = 2^ and ByJp/2 = 
\/l ~ iv. Obviously, it exists for v < 1/3 only. The 
phases <y9i,2 of the fields ■!/'i,2 ~ exp(i(^i_2) change from 
ipi,2 = at a; = —00 to ipi = —ip2 = n a,t x = +co. 

Eq. @ has been studied in a varietyof cornpletely dif- 
ferent contexts and systems jH HI [H EHil • In 
terms of the real and imaginary parts, it was discussed in 
ref. 22] with respect to the transformation of the domain 

walls in magnetics. In the systems mm Hi 113, the 

spontaneous transition between the solutions of eq.@ 
has been found. Its interpretation is system specific. 

DS-.JV interconversion. The DS formally exists for all 
values of the dimensionless coupling v. The JV solution 
is valid only for v < Vc = 1/3. At the critical value Vc, 
the JV turns into the DS. Simple energy argument shows 
that the DS is an unstable state for i' < i^c- The energies 
Eds = f(l + i/)^/2 and Ejy = oftheDS and 

the JV, respectively, as well as their j/-derivatives become 
equal at v — Vc- For v < Vc, one finds Eos > Fjy, 
which implies absolute instability of the DS. 

Despite being identical to the problems of refs. [2^E3 . 
Hi HI m 113 in the static hmit, dynamical eas. p7|) 
cannot be mapped on these systems. Dynamics in our 
case is essentially two-component. The interconversion 
can be well described within the family (jS]), where p — 
p{t), B — Bit) are some real and complex functions of 
time, respectively. Substituting into eas. (|ll5(l and 
performing the variational analysis based on the adia- 
batic approximation with respect to the mass flow along 
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the waveguides, we find for — s- z/c 

d-^b- ^{a^ +b^)b-lb^O, (10) 

i} + - iyr)a + ^{a^ + b^)a = 0, (11) 
v3 

where B ~ a + ib, and a, b are real, and f = 3tj/ (32(1 + 
i^)). It is quite obvious that the interconversion of the DS 
(a = 6 = 0) and the JV (& = 0, a = ±?>^I^^Jv^ - v) pro- 
ceeds on typical relaxation time r of the B J J |^ . Very 
close to the instability ( v ^ the "critical" slowing 
down ~ f l\v — Vc\ takes place. Thus, the DS may van- 
ish in accordance with the mechanism of ref . [3] before it 
decays into the JV. It is important to realize that the 
above result is independent of a particular mechanism of 
dissipation. 

We have also performed direct numerical simulations 
of the full GP equations Ht)l7|) with the initial conditions 
taken as either DS or JV (located at a; = 0) for peri- 
odic boundary conditions, with the space period being 
about 10 times larger than soliton size. To accommo- 
date the phases variation by tt, two-soliton solutions were 
considered. On Fig|21 the results of slow evolution of 
the coupling v from below critical v = 1/7 (where JV is 
stable) toward above critical v — 2/h (where the DS is 
stable) as well as its reverse is presented for the dissipa- 
tion a = 0.5. As perturbation, small uniform imbalance 
of the waveguides population has been imposed on the 
initial conditions. Starting from JV and increasing v 
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FIG. 2: Interconversion of the DS and JV is displayed via the 
evolutions of single waveguide density p{x,t) (solid lines) and 
current J{x, t) (dashed lines) at a; = 0. Slow change of v pass- 
ing through Vc is shown on the inset. Transformation JV^DS 
(DS— >JV) is manifested in the damped oscillatory transition 
of p(0, t) and J(0, t) from finite (zero) to zero (finite) values. 



y, the sign of the JV vorticity is determined by sign of 
the initial density imbalance. We also note that the vari- 
ational ansatz of IjlOl lll|l has been verified to reproduce 
the full numerical solution with good accuracy for the 
initial imbalances < 10%. 

Vortex dynamics. The moving DS solution is known 
analytically (as gray soliton ) [l|. Unfortunately, find- 
ing analytical solution of the moving JV seems impos- 
sible. However, in the limit v the JV can be 
well approximated by the SG equation. Indeed, in this 
case, the variation of the total density can be ignored. 
Thus, the representation ?/'i,2 = v^l ± p/2e^*'^/^ can 
be employed, with \p\ « 1. Substituting this into 
the (dimensionless) eqs.l^-EJ, and ignoring the gradi- 
ent Vp, we obtain after the variation — p — p ~ 
and — p + V^(/3 — AvsuvLp — 0. Without dissipation 
((T — > oo), we find p — if, and, then, the SG equation 
— (/3 + V^(^ — Ai/siiKf — (see in _5j). In the case tr 0, 
one finds p = a(p, and, then, the overdamped SG equa- 
tion [i| —a(p + V'^if — Avsinif = 0. We note that, as 
I' — > 0, the JV solution of eq.JHJ satisfies the static SG 
equation. 

It is important to note that, in the SG approximation, 
the JV is always stable. There is, simply, no room for 
the DS due to the imposed constraints. Similar situation 
has been discussed in ref.||2Sj| for two-component BEC. 

Quasi- ID Fluctuations. At finite temperatures T, 
phase-slip effects can destroy supercurrents. The cor- 
responding life-time, however, can be very long |29l| . In 
our case, stability of the JV is determined by the finite 
energy barrier AE — E^s — Ejy with respect to ther- 
mal (quantum) jumps between the opposite orientations 
of the current circulation. The probability of the thermal 
jump is P ~ exp(— Ai?/r). If — > Vci we find 

P~exp(-4.5V3(i/,-z/)2|l). (12) 

Thus, to have long lived JV, temperature must satisfy 
the condition 



T < 4.5\/3(i/c - vfE* 



(13) 



Here E* = = VpT'c, with = h^nl/m being 
the temperature of the quasi-BEC [2I| formation. For 
J/ — > 0, vortex-antivortex pairs can be created thermally 
(quantum mechanically). To suppress the thermal effect, 
T must be less than the energy 2Ejv of the pair. Thus, 



T < IQ^j^E* 



0. 



(14) 



causes damped phase-slip oscillations during which the 
JV changes its vorticity (see curves marked JV — > DS), 
and, finally, current vanishes. Correspondingly, the den- 
sity, first, acquires zero at finite time moments before the 
zero becomes permanent. This final stage indicates for- 
mation of the DS. If starting from the DS and decreasing 



It is important to note that, if the conditions (|13I14|I are 
met, the JV size Ljy ~ Ic/y^ is smaller than a typical 
phase coherence size L^f, — h^rif)/ [mT] — IcE* /T of the 
quasi-BEC ^30]. 

At low temperatures, quantum tunneling between two 
circulations of the JV will restore the symmetry for 
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1/ < i^c due to the quantum phase-slip. This results in 
a soliton, which, on one hand, is an essentially quantum 
object with respect to its internal structure — a superpo- 
sition of opposite vorticities, and, on the other hand, can 
propagate like a heavy classical particle {Sq ^ 1) along 
the BJJ-waveguides. It can be viewed as a mobile qubit. 

Vortex pump. The vortex can be used to transfer a 
portion of the BEC atoms between BEC reservoirs (see 
FigQ). The rate of the atom deposition/depletion in the 
first reservoir is iVi = Ji, where Ji = pi'^xfi is the cur- 
rent along the first waveguide taken at the point where 
the coupling is zero (the right end of the waveguides). 
It is given by the density pi and the phase <pi at the 
boundary. Thus, AN = dtJi. If the reservoirs are 
large enough, the JV solution can be considered as un- 
deformed at the boundary. Hence, one can set Ji — 
Ji{x — xo{t)), where xo{t) is the position of the JV, and 
replace the integration J dtJi — J{dx/V{x))Ji, with 
V{x) — dxo{t)/dt. If the JV moves uniformly, this gives 
AN — J piS/ifidx/V. For small V, we substitute the 
static JV solution. Thus, the explicit integration gives 
(in the physical units) AiV = ^^%ii^(l -I- i/)(l - 3;/), 
where Vs = fi/jjnlc) is the speed of sound. The GP 
regime 5*0 ^ 1 '6] implies AiV ^ 1. 

Creation and detection of the JV. Observation of the 
described interconversion can be based on, first, creating 
the DS simultaneously in the both waveguides. Then, 
moving them slowly apart (to have v < ly^) will result in 
vanishing of the DS into the JV, so that the zero in the 
densities will heal (see Figl^J. Bringing the waveguides 
back together (to have v > I'c) will cause reappearing 
of the DS (see Figl^J. The very fact of this reversible 
transformation, besides being potentially utilized for cre- 
ating the JV, can serve as an unambiguous evidence of 
the static Josephson currents. It is important that this 
effect, in contrast to the suggestion of ref.i^, can be ob- 
served in the overdamped regime. Direct imaging of the 
JV currents could be done by the Bragg spectroscopy 
technique 31]. Analysis of the absorption imaging of the 
JV upon expansion (for 3D vortices, see [s^) will be pre- 
sented elsewhere. 

Applications. The utilizations of the JV for coherent 
BEC transfer and as a mobile qubit creates quite intrigu- 
ing perspectives for quantum computations and coherent 
BEC manipulations in microtraps. 
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